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Abstract 

Gravitational waves affect the observed direction of light from distant sources. 
At telescopes, this change in direction appears as periodic variations in the apparent 
positions of these sources on the sky; that is, as proper motion. A wave of a given 
phase, traveling in a given direction, produces a characteristic pattern of proper mo¬ 
tions over the sky. Comparison of observed proper motions with this pattern serves 
to test for the presence of gravitational waves. A stochastic background of waves 
induces apparent proper motions with specific statistical properties, and so, may 
also be sought. In this paper we consider the effects of a cosmological background of 
gravitational radiation on astrometric observations. We derive an equation for the 
time delay measured by two antennae observing the same source in an Einstein-de 
Sitter spacetime containing gravitational radiation. We also show how to obtain 
similar expressions for curved Friedmann-Robertson-Walker spacetimes. 

Keywords: Cosmology: Gravitational Radiation, Gravitation- Techniques: Inter¬ 
ferometric 


1. Introduction 

It is commonly agreed that gravitational waves, predicted by Einstein’s theory of 
general relativity (Einstein 1916), must exist in our Universe. To date, however, the 
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only evidence for their presence is that the orbital decay rates of certain binary pnl- 
sars appear consistent with the predicted rates of energy loss from gravitational ra¬ 
diation (Taylor 1992; Taylor and Weisberg 1989). Nevertheless, gravitational waves 
arise generically after the inflationary phase in inflationary cosmologies (Rnbakov, 
Sazhin, and Veryaskin 1982; Fabbri and Pollock 1983; Abbott and Wise 1984) and 
shonld be prodnced in a wide range of physical sitnations at later times (Thorne 
1987; Carr 1980). For these reasons astrophysicists are certain that gravitational 
wave astronomy, thongh difhcnlt, will be of enormons valne in nnderstanding onr 
Universe. 

The effects of gravitational radiation may be divided into two categories. Direct 
effects physically conple the energy density in the waves to matter, cansing, for 
instance, a bar to resonate. Gravitational waves also affect the propagation of 
radiation, cansing a spacetime containing gravitational waves to look different from 
one withont. The very long wavelength (A > 10“^ pc) gravitational radiation 
which we focns on in this paper is best searched for by examining its effects on the 
radiation we receive from astrophysical sonrces. Snch indirect effects have been nsed 
snccessfnlly to constrain the fraction of the energy density in onr Universe which 
can be contained in gravitational radiation of varions wavelengths. Among these 
constraints are Oqw < 10“^ at A ~ 1 pc, from pnlsar timing (Romani and Taylor 
1983; Taylor 1987); Dqw < 0.04 at 10 pc < A < 10 kpc, also from pnlsar timing 
(Taylor and Weisberg 1989); Oqw < 10“^ at A <0.1 kpc if the waves existed dnring 
nncleosynthesis, from nncleosynthesis constraints (Carr 1980); Oqw < 10“® or 10“^ 
at A > 1 Mpc if the waves did or did not exist, respectively, at recombination, from 
microwave backgronnd anisotropy limits (Linder 1988a); and Oqw < 10“^ for 30 
kpc < A < 300 Mpc, from galaxy-galaxy n-point correlation fnnctions (Linder 
1988b). 

Recently Enbanks and Matsakis (1994) have reported Very Long Baseline In¬ 
terferometry (VLBI) measnrements that indicate qnasars have a dehnite pattern of 
apparent motions on the sky with root-mean-sqnare (RMS) angnlar velocity ~ 20 
/uas yr“^. The work of Linder (1988b) fnrnishes an estimate of the constraint we can 
expect from data of this accnracy. Linder obtains an expression for the mean-sqnare 
angnlar deflection of light from cosmological sonrces indnced by gravitational waves. 
Dividing this expression by the sqnare of the wave period we gain an estimate for 
the mean-sqnare angnlar velocity of sonrces at redshift ^ in an Einstein-de Sitter 
spacetime indnced by waves with energy density as a fraction of the closnre density 
Oqw given by ~ (^1 -f (1 -t- H^Qgw where Ho is the Hnbble constant at 

z — 0. Using Ho = lOO/i km s~^ Mpc~^ this tells ns (taking ^ = 1 as a representa¬ 
tive redshift) that the qnasar motion data may be expected to either constrain or 
detect gravitational waves at a level of Oqw = 0.04h“^. This is competitive with 
the pnlsar timing limits, and shonld cover a mnch larger range of gravitational wave 
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wavelengths. 

In order to test the hypothesis that the qnasar motions reported by Enbanks 
and Matsakis (1994) are cansed by a cosmologically signihcant backgronnd of grav¬ 
itational waves we need a theoretical framework snitable for the analysis of the 
effects of snch waves on VLBI measnrements of distant sonrces. In a seminal work 
on observations in cosmology, Kristian and Sachs (1965) established a nnmber of 
formnlae relating the observed properties of cosmological sonrces to the physical 
properties of the spacetime in which they are observed. Their formnla for the 
proper motion distance conld be nsed immediately to analyze the system we are 
concerned with here except that their work ntilizes an expansion in the distance to 
the sonrce divided by some reasonably dehned radins of cnrvatnre of the spacetime. 
For the high redshift qnasars of Enbanks and Matsakis (1994) snch an expansion is 
not nsefnl as the qnasars are a large fraction of the Hnbble distance, and many grav¬ 
itational wave wavelengths, from ns. The work of Linder (1988b) also bears close 
relationship to the problem nnder consideration. Linder has obtained the deviation 
in angle snffered by a light ray in an Einstein-de Sitter spacetime containing grav¬ 
itational waves. For astrophysical thin lens systems, where the angnlar deflection 
of an incident light ray may be considered to occnr at a single point, knowledge of 
the angnlar deviation is snfhcient to determine the apparent position of the sonrce 
on the observer’s sky. This is the content of the well known lens eqnation. The 
sitnation is different, however, when a cosmological backgronnd of gravitational 
waves is effectively acting as the lens. For the wave case, the angnlar deflection 
occnrs over the entire photon path and there is no obvions a priori relationship 
between the sonrce position on the observer’s sky and the pnrely mathematical, i.e. 
nnobservable, angnlar deviation. 

In this paper we nse the pertnrbative geodesic expansion introdnced in Pyne 
and Birkinshaw (1993) to determine the effects of gravitational radiation on VLBI 
measnrements of distant sonrces. In Gwinn et al. (1995) we nse the resnlts of this 
work to test the hypothesis that the qnasar motions reported by Enbanks and Mat¬ 
sakis (1994) are cansed by a cosmologically signihcant backgronnd of gravitational 
waves. 

The ontline of this paper is as follows. In section 2 we develop a method for 
analyzing a VLBI experiment in metric pertnrbed Einstein de-Sitter spacetimes. 
In section 3 we apply the method to the case of an Einstein-de Sitter spacetime 
pertnrbed by a spectrnm of cosmological gravitational waves. In section 4 we con¬ 
sider a single, plns-polarized, monochromatic wave and determine the pattern of 
sonrce proper motions which it prodnces on the sky. In section 5 we generalize onr 
eqnations to the cnrved Friedmann-Robertson-Walker (FRW) spacetimes and show 
how a physical nnderstanding of the method emerges from an analysis of the Jacobi 
eqnation. In section 6 we investigate the consistency of onr eqnations by considering 
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a simple gravitational lens system. In section 7 we present onr conclnsions. 


2. Light Rays in a Perturbed Einstein-de Sitter Spacetime 

The resnlts of Pyne and Birkinshaw (1995) allow constrnction of the paths 
of light rays throngh a pertnrbed Einstein-de Sitter spacetime with only minimal 
effort. The metric for snch a spacetime takes the form 


+ dy^ -t- dz^^ -|- a^h^i.dx^dx’^ (1) 


where a is the Friedmann expansion factor. We let = o^ds^. By standard 
conformal resnlts, light rays in ds^ and ds^ coincide and their (affine) parameteri- 
zations are related by = a~‘^k^. Here, and thronghont this paper, Roman letters 
ijj,--- rnn over {1,2,3}, while Greek letters /U, z/,... rnn over (0,1,2,3}. We nse 
geometrized nnits, G = c = 1. The spacetime metric is taken to have signatnre +2. 
Onr Riemann and Ricci tensor conventions are given by [Aa, 2S.p]v^ = 
and Rdf) R^dfifd- 

The metric, (1), is of a class of metrics whose radial nnll geodesics were investi¬ 
gated in Pyne and Birkinshaw (1995). For the specihc case of (1), the Einstein-de 
Sitter backgronnd, non-radial nnll geodesics may be constrncted from the resnlts 
of that paper with almost no effort (this is becanse the Jacobi and parallel prop¬ 
agators for these geodesics are the same as those nsed in that work. For the case 
of the cnrved FRW backgronnds, the propagators for the radial and the non-radial 
geodesics are not eqnivalent). We express this in the form of a 

Theorem: Let denote that part of ds^ independent of the pertnrbation (i.e. 

the Minkowski metric). Let be an affinely parametrized nnll geodesic of 

with = dx^^^^{X)/dX snch that k^^^^ = 1. Pnt 


/(i)m = -r(i)M^^^(o)a^(o)/3_ 


( 2 ) 


with that part of the Christoffel co nn ection of ds^ which is linear in the 

metric pertnrbation and its hrst partial derivatives. Then x^{X) — x^^'^^{X) -|- 
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ig an affinely parametrized geodesic (not necessarily nnll) of ds^ to first 
order provided that, for all A 2 , Ai 


X 


(i)m 


(A2) = (Ai) + (A2 - Ai) (Ai) + (A2 - A) /(')^(A) dX ( 3 ) 

J\\ 


where = dx'^^^^{X)/dX and the integration is performed over a:*^°^^(A). 

The condition = 1 on the affine parametrization of the backgronnd geodesic 
in the above theorem is imposed simply so that we can nse the propagators of Pyne 
and Birkinshaw (1995), who imposed that condition for ease of calcnlation. It is not 
hard to compnte the necessary propagators for with r any affine parameter. 

This is not really needed, however, since the geodesic, k^(X), constrncted by the 
theorem above may be reparametrized directly. For this reason we will sometimes 
refer to k^{X) as a wavevector thongh this term is nsnally reserved for the tangent 
to a nnll geodesic parametrized so that u^k^ is the photon freqnency observed by 
an observer with fonr-velocity u^. 

Snppose now that we solve (3) along some given geodesic of x^'^^^{X) 

snbject to x^^^^ (A 2 ) = 0 and x^^^^ (Ai) = 0. We can not simply demand that the 
separation, a:^^^^(A), vanish at both Ai and A 2 if we want the constrncted geodesic to 
be nnll since x^^^^ (Ai) and x^^^^ (A 2 ) are nnll separated in bnt not necessarily 

in ds^. Becanse the spatial and timelike components of (3) deconple, however, we 
can nse the above bonndary conditions to obtain (Ai). We can then solve for 
(Ai) by demanding that onr constrncted geodesic be nnll. 

The condition that onr constrncted geodesic be nnll in can be written 


^( 1)0 _ ^ 
2 


( 4 ) 


where the eqnation holds along (A). Taking onr bonndary conditions into 
acconnt, this allows ns to write the timelike component of (3) as 
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a:WO(Ai) 


( 5 ) 


= -1 (A2 - Ai) - (A2 - Ai) 

- r'(A2-A)/W0(A)dA, 

J Ai 


the inner prodncts being evalnated at a:^°^^(Ai). The spatial components of (3) 
yield 


(A 2 - Ai) (Ai) - - (A 2 - A) /(i)*(A) d\ 

JXi 


( 6 ) 


which may be combined with (5) to prodnce 


(Ai) = (A 2 - Ai) (A 2 - A) dX (7) 


Eqnation (7) conld also have been obtained immediately from (3) and (4) after 
taking the inner prodnct of (3) with 

At this point another representation of the pertnrbation vector, is very 

nsefnl. Letting a semicolon denote covariant differentiation nsing the Christoffel 
connection of we have 


/(1)m = ^/^^^;M^(0)a^(0)/3 _ 


( 8 ) 


Since is geodesic, this gives 
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(9) 


This may be substituted into (7) allowing an integration by parts to produce 


(Ai) = -i dX 


( 10 ) 


We have thus found the following 

Corollary: Given two points, q = {X 2 ) and w = (Ai), connected by 

a null geodesic of with = 1, then to hrst order the points 

w' = (Ai) + (Ai), 0*^ and q are null separated in ds^ along some geodesic 

x^ provided (Ai) obeys (10). Further, the tangent vector to x^ at w' is set by 
(4) and (6) above. 

In Section 3 we will show how this corollary may be used to analyze the effects of 
an arbitrary metric perturbation on VLBI observations in the metric (1). 

We hasten to point out that the theorem, and so corollary, above have no 
more content than a direct integration of the linearized geodesic equations of our 
spacetime. We have presented this information in this ma nn er for two reasons. 
First, as we will demonstrate in the next section, the above corollary is specihcally 
adapted to an analysis of VLBI experiments in the spacetimes we are considering. 
Second, the presentation above is organized so as to facilitate the generalization 
to curved backgrounds described in section 5 below. The reader will note that for 
scalar perturbations in the longitudinal gauge (see e.g. Mukhanov, Feldman, and 
Brandenberger 1992), the RHS of (10) is the Shapiro delay evaluated along the 
background path. We will return to the relationship between (10) and lens systems 
in section 6. 

It remains for us to discuss the consistency criteria for the manipulations leading 
to our corollary. We defer an examination of this topic till section 6, contenting 
ourselves here with an informal remark. Imagine a geodesic of which 

intersects w' at affine parameter value Ai with tangent vector + k^^^, 

where a subscript w {w') denotes evaluation at w {w') and is set by (4) and 
(6). We can use this geodesic in (3) to construct a null geodesic in ds^, x^{p), 
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obeying 5;^(Ai) = w' and At a minimnm, then, we wonld expect 

(10) to be a consistent solntion of the hxed endpoint problem only if (A 2 ) = q. 
Ronghly speaking, we wonld expect this to hold if the gravitational effects of the 
pertnrbation are similar on and for eqnal affine parameter valnes, that 

is if X and x pass throngh “snfhciently similar” metrics. 


3. Gravitational Radiation and VLBI 

In order to nse (10) to investigate the effects of a gravitational wave backgronnd 
on a VLBI experiment we hrst need the appropriate form for the metric pertnrba¬ 
tion. The plane, monochromatic wave solntion in the synchronons gange to the 
pertnrbed Einstein’s eqnations associated with (1) can be written as the real part 
of 


hoo 

h-Oi 

hij 


0 

0 


h+ (p) ^Re+R^'^,, + hy^{p) [Re^R^^,, 

IJ 


^i{p-x-pri) 


( 11 ) 


(Hawking, 1966), where ag is some hdncial valne of a, h-t- and hx are complex 
valned fnnctions containing the amplitnde and phase information, e_(_ and Cx are 
polarization matrices which we represent by 


e+ 


ex 



0 0 \ 
-1 0 
0 oj 

1 0 \ ’ 
0 0 
0 oJ 


(12) 


J? is a rotation matrix, R^ its transpose, which we represent by 



R{OA) 


sin (f) cos 9 cos (f) sin 9 cos (f) 

— cos (f) cos 9 sin (p sin 9 sin (p 

0 — sin 9 cos 9 


(13) 


and the modevector, 


( sin 9 cos (p \ 
sin 6* sin 
cos6' ) 



(14) 


We can think of the angles, 9 and (p, as fnnctions of p throngh (14). This solntion is 
valid in what is known as the adiabatic regime, which demands that, at the times 
of interest, the physical rednced wavelength of the wave is mnch smaller than the 
Hnbble distance. Mathematically this is simply 1 ^ a'/p where ' = d/dt, t the 
comoving time co-ordinate related to r] by d/dt = a~^d/dr]. The comoving rednced 
wavelength, 1/p, is related to the physical wavelength on a spatial hypersnrface 
of constant conformal time ry, Aphys, by Aphys = 27ra{r])/p. The condition that 
a wave be in the adiabatic regime at (nnpertnrbed) redshift z is thns Aphys ^ 
27r(l -|- , with Ho the Hnbble constant at z = 0. The general tensor 

pertnrbation in the adiabatic regime, is a snperposition of the plane, 

monochromatic waves for different modevectors. 

We now show how to nse the corollary of section 2 to analyze a VLBI exper¬ 
iment in a metric pertnrbed Einstein-de Sitter spacetime with pertnrbation given 
by (11). The basic idea is to constrnct the fntnre nnll cone of a point of emission by 
determining the time at which the cone intersects any given line of constant spatial 
position. This will allow ns easily to determine the difference in time of reception 
of a signal at two points a known spatial distance apart, that is, at two ends of an 
interferometer. The geometry of this section is shown in Fignre 1. 

Snppose an antenna, which we label antenna A, receives at the point w'j^ = 
(VaHa) photons from a sonrce at q = {rjo — L^Lcs). We take Cg to be normalized 
to nnity in ds^^R. Denote the translation of w'j^ throngh time by ta = {fj, ta)- The 
forward nnll cone in ds^^R of q intersects ta at the point wa = (w + ta) where 
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Aa = \Les - taI - L. 


(15) 


In (15) and below, the inner prodnct is written nsing vertical bars. We 

will also nse a raised dot for this inner prodnct. We note that Aa has a very 
simple interpretation: two antennae, one at the spatial origin and one at spatial co¬ 
ordinates given by ta, observing a sonrce at q in an nnpertnrbed Einstein-de Sitter 
spacetime wonld measnre a time delay of Aa between the reception of a signal at 
the origin at time rjo and the reception of the same signal by the antenna at ta- 
The connecting nnll geodesic is given by a:^^^(A) = (A, Leg — {X — rjo + L) ca) with 

Leg - Ta 

CA = jy-. -^ 

\Leg - rA 



which is seen to be properly normalized. The corollary of section 2 now demands 


7j'a = 7Jo + Aa- 1 r dX (17) 

2 JtIo+Aa ^ ^ 

the integral being taken over a:A^^(A). 

It will be convenient for ns to nse not (17), bnt the eqnivalent expression in 
terms of the comoving time. Provided that the change in the scale factor over Aa 
and ? 7 a^ may be neglected 


f A — fo + ^qAa + ^A^ 
ttn 


= to + ttoA 


A 


2 J77o + Aa 


dX, 


the integration occnring over a:^^^(A). 


(18) 
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We now consider a second antenna, antenna B, which receives photons at = 
from the same point of emission, q. This system obeys eqnations (15) and 
(18) when all of the snbscript A’s are replaced by snbscript S’s. The measnred 
(comoving) time delay for the two antenna system is given by 


Td — t's 

= tto {Aa - Ab) + t^A - tB^ 


(19) 


with tA^ given by the term in (18) linear in h^i, and given by the analogons 
expression for antenna B. 

To see that (19) contains the nsnal expression for the time delay of a VLBI 
system, consider the term, in (19) which is independent of the pertnrbation. 

Pnt tb = TA + b = r + b. We take the distance to the sonrce to be mnch larger than 
the proper lengths of either r or 6, and we work to hrst order in the implied small 
qnantities, e.g b/L. This is the mathematical expression of the locally plane wave 
approximation whereby we neglect the cnrvatnre of the wavefronts at the observer. 
With this approximation we have ~ ca and 


do (AA — Ab) 


^ Gob ■ Cs 


( 20 ) 


Continning, we note that the scale factor Go ^ a {t'A) = ci(t) where we have written 
the time of reception at A, I'a, simply as t. Provided antenna B does not move sig- 
nihcantly over Td we can write b as b{t ), the (comoving) co-ordinate vector reaching 
from antenna A to antenna B at time t. In addition, we can give meaning to es{t) 
as follows: the spatial co-ordinates of the sonrce at the event of emission of the 
photons which arrive at antenna A at time t is given by L{t)es{t), with es{t) nor¬ 
malized to nnity in The hrst term in (19) can then be written gE ■ Cg, with 

each qnantity evalnated at t. 

The tensor which projects the wavevector of the arriving photons into the 
rest space of the observing antennae is written 6^ + with the antennae 

fonr-velocity. Taking onr antennae comoving (for simplicity, a Lorentz boost will 
easily yield the general case from the specihc) the properly normalized (in ds^) 
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four-velocity is given by — (l/a,0). This writes the unit (in ds"^) vector in the 
rest space of the observer which points toward the source, as 


s 


M — 


u~k 



( 21 ) 


It is easy to check that has no timelike component. In (21) we have used ■ to 
denote the inner product. Since k^ = a~‘^k^ = -|- a we have 


as 


* = e* - 


— e 




( 22 ) 


where we have used (4) and we have written ca simply as e. At zeroth order this 
allows us to write (19) as 


rp(0) 



(23) 


This formula, scalar under spatial transformations in the antenna rest space and 
whose quantities are simply the restriction to this space of tensor quantities (which 
by dehnition have well dehned properties under Lorentz transformation), is exactly 
the standard formula of VLBI (Thompson, Moran, and Swenson 1986). 

Under the locally plane wave approximation, (b'^g^jS^^ is the spatial distance in 
the observer’s rest frame that the incoming wavefront must still travel to antenna B 
at the instant it has hit antenna A. It is, thus, the time delay. In this approximation, 
then, we can obtain the time delay either by determining the perturbed direction 
vector to the emitter, s*, and expanding its inner product with the baseline vector, 
6*, or we can work directly from (19), imposing the mathematical constraints which 
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lead to locally plane waves at the observer. We now show that these methods agree 
by calcnlating the explicit formnla for the time delay nsing both methods 

We have already seen that the zeroth order expressions agree. Unfortnnately, 
the proof to hrst order is not so simple. We start by setting ao = 1 in order to 
simplify onr notation. This means that at wa there is no distinction between ds'^ 
and ds'^. Since we will only be concerned with qnantities at wa in onr formnlae, 
we will, for the rest of this section, no longer pnt a bar over the physical metric. 
Nevertheless, its presence is implied, that is, e.g. may be nnderstood as 
the two being eqnal at the point we need them. 

We hrst calcnlate the hrst order term in the time delay from the formnla 


T^(l) _ 

— 



(24) 


From (22) we see that we need the pertnrbation to the spatial components of the 
wavevector at antenna A, k^A^. This qnantity is solved for by (6), above, nsing the 
backgronnd path with the resnlt 


k 


{i)i 

A 


. 2 

Qca ■ ddeAP -ip(l+p- ca) Hca^ ^ '' 

^ ^ip-eAiVo+AA) J 


(25) 


with 


J = 


(Po-L) r ^ r ^ dX (26) 

Jrjo+AA Jv^+Aa 


Vo-L 

T/o+Aa 


The integrals may be performed explicitly, resnlting in combinations of exponen¬ 
tial and exponential-integral fnnctions of argnments —ip{l +p ■ ca) (vo ~ L) and 
—ip {1 + p ■ ca) ivo + ^a)- The resnlt, while correct, is not particnlarly easy to 
nse and we do not write it down. Instead, we note that adiabaticity demands 
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p{'r]o — L) 3> 1, which allows us to use the large argument expansion for the 
exponential-integral functions (e.g. Gradshteyn and Ryzhik 1994, equation 8.215) 
for all values of ca not too close to —p. Next we suppose that our source is many 
gravitational wave (reduced) wavelengths away, pL ^ 1, and we agree to work only 
to leading order in {j>L)~^. The solution, for the appropriate range of e^, then 
becomes 


.■ = 8!,+ , 

(?7o A a) 






(27) 


Because this expression shares with the exact result, (25) and (26) above, the prop¬ 
erty of vanishing at ca = —p we will use it for all values of ca- Rigorously, we 
should examine its approach to this zero and compare it with that of the exact 
expression, but we will consider ourselves justihed by the end result. Noting that, 
to hrst order. 


hij {w'a) = hij {wa) 

2 

_ Vo _Qip-rAQ-ip{rtO+^A)pp 

(tjo + A^)^ 


(28) 


equation (27) writes for the time delay 


Td = 


= h-eA + -e\h 


1-J^A 


b ■ eA + b- p 
{I+P-ca) 


(29) 


We note that this formula is a well-dehned scalar in rest space of our observer. 
We also point out that p = Pphys and that our hatted vectors are of unit norm 
in the physical background metric at the observer since, briefly returning to our 
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old bar notation, = rj^^. Defining the dimensionless vector ^ + p, with 

C = C/^C^^Cb we have 


Ta = h■eA + h■C(Ch^Jl^). 


(30) 


In the following section we will show how this formnla may be nsed to obtain the 
pattern of proper motions on the sky measnred by an interferometer observing a 
sample of distant sonrces throngh a backgronnd of gravitational radiation. 

For completeness, we now show how onr solntion, (30), is obtained from (18) 
and (19), above. These write the first order contribntion to the time delay as 




t 


(1) hi) 


A 


B 


1 fVo-L 


2 JtIo+Aa 






{x^a) k 


r<JA + i 


rVo-L 

Irjo+AB 



dX. 

(31) 


In this formnla, the path on which hfj,,y is to be evalnated has been shown explicitly. 

The easiest way to proceed is to expand the integrand in the integral for 
abont the valne of the integrand in the integral for for fixed A. If we neglect 
terms of qnadratic order in b/L, 


k 


(0)m 

B 



k^A^h^^ (x^a) k^A'" + 25k^^^^h^^ (x^a) k^A 
+ k^A^^h^jy^a (x^a) 


(32) 


with 


= k^B^ - k^A^ 

= -( ^ _p \ (33) 

L \b- (6-e,)eJ 
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and = {X + L — rjo) Using (32) in (31) yields 


= 1 r ^ V (x^a^) dX 

2 Jr,o+AA ^ ^ ^ ^ 

+ V (pa) k^A^' 

= \kA^^h^^ {pa) k^A'" (34) 

2 

-eA-H(b-{b- es)es) 

^ ' Jj 

2 

- ■ HcaP- {b- (6 ■ e,)e,) ^e^v■rA^-^v■eA{m+AA)J 


where J is given by (26) and 


/= r“ ^ A-2e-*P(i+P-^^)^dA. 

Jrjo+AA 


Once again, the integrals we need to perform yield a combination of exponen¬ 
tial and exponential-integral fnnctions of argnments —ip{l +p ■ ca) {rjo ~ L) and 
—ip (1 -|- p ■ ca) {rjo + A^). We follow onr earlier conrse and nse the large argnment 
expansion of the exponential-integral fnnctions. The terms in the resnltant expres¬ 
sion contribnted by the term written proportional to / in (34), are of order 1/pL 
times the terms contribnted by the term written proportional to J in (34), and of 
order 1/pL times the hrst term in (34) (in fact, direct dimensional analysis of (34) 
argnes for this conclnsion). Again restricting onrselves to sonrces which are many 
gravitational wave rednced wavelengths away, we may neglect the terms down by 
pL. The remaining terms combine easily to give the hrst order part of (29), which 
is what we needed to show. That is, we have proven, throngh hrst order. 


tA-ts ^ b’^QijS^ 


(36) 
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with the LHS assembled from (15), (18), and (19), and the RHS constrncted from 

( 22 ), 

Becanse this has been a long section we will collect the approximations nsed 
to gain onr formnla here. We envision an adiabatic backgronnd of gravitational 
radiation and two antennae, A and B, snch that the co-ordinate separation vector, 
6, reaching from antenna A to antenna B has a proper length mnch smaller than the 
(rednced) wavelengths of the waves, pb 1. We snppose the antennae to be nearly 
comoving over the timescale of the measnred delay. We farther snppose that the 
sonrces observed are many gravitational wave (rednced) wavelengths away, pL 3> 1. 
Provided that these ineqnalities hold, the measnred time delay, is given by (29), 
with ca related to the observed direction vector to the sonrce, s*, by (22) (nsing 
(25) as well). We note that all of onr conditions may be relaxed if necessary, by 
nsing the more precise eqnations (19) and (34). 

We can obtain a very rongh a posteriori consistency check on the accnracy of 
onr pertnrbative expansion in the following way. We approximate the separation 
between the solntion geodesic, (r^(A), and by (Ai) ~ Lh, where h is 

some characteristic element of That is, the separation is ronghly the distance 
to the sonrce times the angle between the two geodesics at the observer. We mnst 
demand that the pertnrbations “felt” by these two geodesics are nearly identical, 
which means, since varies over scales larger than 1/p, Lh -C p~^. The idea is 
now to sqnare this ineqnality and relate to the fraction of the closnre density 
contained in the gravitational radiation. An explicit model for the wave backgronnd 
is necessary for this step. For illnstration, consider a stochastic backgronnd of 
gravitational radiation. Then Oqw ~ p^h^/i^o- Since 


L = {l - (1 + (37) 

we can write onr constraint as 

k^GW < ^ (l - (1 + ^) (38) 


This is to be interpreted as follows. If we nse the formnlae above to infer, from 
sonrces at characteristic redshift z, the presence of a stochastic backgronnd of grav¬ 
itational waves with an energy density of flow, then a crnde estimate of consistency 
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is furnished by the degree which (38) is found to hold. We point out that for ^ ~ 1 
this is not a highly restrictive condition. 


4. The Pattern of Proper Motions on the Sky 

In this section we will show how equation (30) may be used to gain the pattern 
of proper motions inferred by an interferometrist observing distant sources through 
a background of gravitational radiation. Since we are working in the linear regime, 
without loss of generality we consider the effect of a single plus-polarized monochro¬ 
matic wave. The effect of a general wave background may be found from the results 
of this section using superposition. We choose our coordinates such that the 2 ;-axis 
is aligned with the waves direction of propagation. The modevector, p, is then given 
hj p = pz and the metric perturbation is written 


hoo 

hoi 

hij 


0 

0 

(e+) 

a 


ij 


(39) 


As in the last section we have set ao = 1. Unlike there, however, in this section 
we return to our practice of using an overbar to distinguish between the perturbed 
Einstein-de Sitter spacetime and its conformal relative, shown in (1). 

The approximations of the previous section guarantee the existence of a region 
of spacetime E such that E contains the events of reception, a ~ 1 in E, and px^ 1 
for any G E. This last condition simply states that the spatial extent of E on 
any constant conformal time hypersurface is much smaller than the wavelength of 
the gravitational wave. In E, to within our level of approximation, the coordinate 
transformation 
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(40) 


t' =t 

i' = X 

</' = (l - ih+e-Pl"-”*) y 


brings the actual metric to Minkowski form; We note that in (40) r] 

is to be considered a function of t in the usual way. 

For simplicity we again suppose at least one of the antennae is comoving, having 
four velocity = (1,0, 0,0). The surfaces of constant conformal time then serve 
as instantaneous rest three-spaces for an observer at this antenna. The co-ordinate 
transformation (40) may be considered a transformation of co-ordinates in these 
hyper surfaces. Since the spatial components of s*, given by (21) above, dehne a 
vector in these hypersurfaces, the observed direction to a source in the primed 
coordinates is simply 



(41) 


and the time delay measured by an antenna pair with separation 6* is given by 



(42) 


Consider now the triad of vectors 
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f-d) = B (l - i'l+e-*'”') i 

»(2) = B (l + ift+e-’”') S 

6(3) = Sz 


with X, y, and z the unprimed (spatial) co-ordinate basis vectors (so, for example, 
= (0,0,1,0) ). The vectors 6 (j) are orthogonal and have proper length B with 
respect to Their spatial components dehne a triad of three-vectors in the 

observers rest space, orthogonal and of proper length B with respect to the induced 
metric Qij. In terms of the primed (spatial) co-ordinate basis vectors, x', y', and z\ 
the triad is written 


6 ( 1 ) = Bx 

b{ 2 ) = By’ (44) 

6 ( 3 ) = Bz' 


We now suppose that the observer makes time delay measurements of a source 
with ca = C(X + Py + jz using each of the three-vectors of the triad as a baseline. 
Then (30), (39), and (43) tell us that T(i(j), the time delay measured along baseline 
6 (i), is given by 


Td{i) 

Td{2) 

Td{3) 


= aB 


^PB 


(l - ) 

(l + Ift+e"''”') 


1 + Ift+e-’’”' 
2 

1 + -h+e-’P” 

2 


= B 


7 + h+e-*'”' (a" - f) 


(Tf)l 


(45) 


Since the time delay is a scalar under co-ordinate transformations in the observer’s 
rest space we can also compute the delays in the primed co-ordinates. 
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(46) 


Tdi^) = h\ 


:)Vr. 


r'S 


Using (44) in (46) and eqnating the resnlting expression with the right-hand side 
of (45) allows ns to read off the components of the normalized direction vector to 
the sonrce in the primed co-ordinates, 


- a (^1 - 
sv' = /3 (l + 


1 -F ]-h+e-^P^ 

1 -f -h+e-^P^ 

2 




7 + 


(a^ - 


m] 

(=rif)l 


(47) 


The proper motion inferred by the observer is given by = P* p/u'^ |q,/ 

where P“ p/ — -|- u°‘ up/ is the projector into the observer’s rest space and a 

slash denotes the covariant derivative of the pertnrbed metric, ds^. To first order 
this resnlts in the simple = s* A comoving sonrce will have a, /3, and 7 
constant. Then 


/2 


ipa 


h-^e~ 

ip(3 

ip 


■ipr) 


a 




X 


h+e-^P^ 1-F 


1 + 7 
o? — (3^ 


1 + 7 

X-h+e-*'”’ (a^ - /3^) z' 


(48) 


Allowing the direction cosines to depend on time enables the effect of sonrce pecnliar 
motion to be calcnlated. To the level of accnracy considered here the effect of sonrce 
motion is to snperpose the standard FRW proper motion resnlts on the gravitational 
wave pattern given by (48) (see e.g. Weinberg 1972, Chapter 15). 
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The expression, (48), for the proper motion contains the direction cosines a, (3, 
and 7 which are the components of the nnobservable vector ca in the nnprimed co¬ 
ordinates. These differ from the direction cosines in the primed co-ordinates, a\ (3', 
and 7 ' where ca = O-'x' + (3'y' , by terms of hrst order. Fnrther, as can be seen 

from (47), a, /3, and 7 describe the components of the observed direction vector to 
the qnasar, s, to zeroth order in the primed (and nnprimed) co-ordinates. Since jl 
is entirely hrst order we may reinterpret a, /? and 7 in (48) as the direction cosines 
to the sonrce in the primed co-ordinates, the difference between these qnantities 
being hrst order and so resnlting only in an ignorable second order correction to 
jl. In terms of the spherical polar co-ordinates of the sonrce in the primed frame, 
dehned by s = sin 9' cos (f)'!:' -|- sin 9’ sin (f)'!/' -\- cos 9'z'^ we have 


a ~ sin 9 ' cos (f)' 

(3 ^ sin 9'sin (f)' (49) 

7 ~ cos 9' 


the approximate sign denoting eqnality to zeroth order. Snbstitntion of (49) into 
(48) and resolntion of the resnltant expression into components along the spherical 
polar basis vectors given by 


9' = cos 9' cos (l)'x' -t- cos 9' sin (f)'!/' — sin 9'z' 
— sin (l)'x' -|- cos (j)'y' 


(60) 


prodnces the resnlt 


jl = sin 9' (cos 2 (/)' 9' — sin 2 (/)' 


(51) 


expressing the inferred proper motion in terms of the wave parameters, p and /i_|_, 
and the angnlar co-ordinates of the sonrce in a gaussian normal frame (to the 
necessary order of approximation) at the observer. 
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5. Generalization to Curved Backgrounds 

In this section, we will consider the general context for the calculations of the 
preceeding sections. By using the Jacobi equation of the background spacetime, 
we will be able to obtain a physical description of the manipulations involved. We 
begin by recalling that the perturbative geodesic expansion (PGE) constructs an 
affinely parametrized geodesic of a general perturbed metric, from an 

affinely parametrized geodesic, of the background metric, by writing 

the sought after geodesic as x^{\) = and solving for 3 ;(i)m(A), the 

separation (Pyne and Birkinshaw 1993). Let the Jacobi propagator along 
be written in terms of its 4 x 4 subblocks as 



^ (Ai, A 2 ) 
C'(Ai,A2) 


D{\iA2))- 


(52) 


Then the equation governing the separation is written 


P (Ai, A 2 ) x^^^ (A 2 ) — A (Ai, A 2 ) x^^^ (Ai) + B (Ai, A 2 ) (Ai) 


+ B (Ai, A 2 ) 




(Ai 


A=Ai 


(53) 


+ /^'s(A2,A)P(Ai,A)/W (A) d\ 

J Ap 


where = dx^^^/d\, P is the parallel propagator along a:^°^^(A), and is the 
perturbation vector given by (8), above. In (53), we have employed a 4 x 4 matrix 
notation so that, for example. 


/>(Ai,A)/<‘'(A)=P(Ai,A)'‘„/<'>"(A) . (54) 


Specihc forms for the Jacobi and parallel propagators of the curved FRW spacetimes 
can be found in Pyne and Birkinshaw (1995). 
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Equation (53) is the generalization of equation (3) of section 2. The proper 
generalization of equation (4) of section 2, the condition that k^{X) be null to hrst 
order in the perturbed metric, is written 


= 0, (55) 


the constraint holding at each point of a:^^^^(A) provided it holds at any given 
point. In principle, equations (53) and (55) allow us to carry out analysis of the 
geodesic problem subject to hxed-endpoint (or mixed) boundary conditions in arbi¬ 
trary metric perturbed spacetimes. Such boundary conditions are more applicable 
to certain astrophysical systems, such as multiple image leasing or VLBI, than the 
more common initial-value boundary data which is useful, for instance, in studies 
of the Sachs-Wolfe effect (Sachs and Wolfe 1967). 

The primary use we will make of equations (53) and (55) in this paper, how¬ 
ever, is to understand the manipulations of sections 2 and 3. The physical picture is 
made clear by recognizing that (53) is exactly the Jacobi equation of the background 
spacetime subject to a forcing perturbation Consider the hxed endpoint so¬ 

lution between q and which made use of the unperturbed geodesic 
between q and wa- We will simply state the conclusion; the reader can carry out 
the steps explicitly and compare them with the development in sections 2 and 3 
to conhrm the conclusion that we offer here. The solution is constructed in the 
following way. First, solve, using (53), for the spatial separation, 5 x{\2)^ on the 
constant conformal time hypersurface containing the point of emission, q, attained 
by a perturbed null geodesic which intersects wa with wavevector coincident with 
/j(o)m (Ai). We can think of this as tracing the photon into the past and determining 
its intersection with the constant conformal time hypersurface containing q. The 
proper pertubation to the spatial components of the wavevector at wa, (Ai), is 
then taken to be that perturbation at wa which, when considered as an impulsive 
perturbation in the background Jacobi equation, produces a deviation vector on 
the constant conformal time hypersurface containing q equal to dx (A 2 ). Speaking 
in informal language borrowed from gravitational lens theory, we solve a standard 
initial-data problem to gain a spatial separation in the source plane, then we con¬ 
vert that to an angle at the observer using the angular diameter distance of the 
background. 

From the spatial components, (Ai), we determine the timelike component 
of the wavevector perturbation, /c*^^)^(Ai), by imposition of the null constraint 
(55). Imagine we now constructed a null geodesic of g^J which intersected wa 
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with wavevector (Ai) (Ai) and used this geodesic to obtain a perturbed 

geodesic which intersects wa with coincident wavevector. Modulo questions of con¬ 
sistency, which we will discuss below, the resultant perturbed geodesic will intersect 
a time translate of q but not, generally, q itself. That is, it will “hit” the proper 
spatial co-ordinates but will “miss” the point q by some offset in time. We correct 
for this by “moving” the point wa in time to another point w'a- The necessary time 
translation is (Ai). 

It is clear that great simplihcations have resulted from the foliation of our back¬ 
ground spacetime by spatial hyper surfaces. While the FRW backgrounds possess 
such a foliation a general background does not, and we have not considered hxed- 
endpoint problems in such spacetimes. In the next section we return to general 
questions of consistency, using as a theoretical laboratory the classic astrophysical 
example of a hxed-endpoint solution, the gravitational lens. 


6. A Simple Gravitational Lens and Consistency 

We work in a Minkowski space background. We consider an observer at the spa¬ 
tial origin of co-ordinates, a lens, of mass m, located at spatial co-ordinates Lx, and 
an emitter which emits a burst of photons at spacetime event q — (—2L, 2L, 0, 0). 
Figure 2 illustrates our geometry. The question we want to answer is where on the 
sky does the observer see the emitter? The answer, of course, is that the observer 
sees a circular ring around the lens, with the angle between lens and ring given 
by the Einstein angle, 6 e = \j2m/L (in our example the ring appears only for an 
instant, but this is unimportant). We want to see how this result emerges from 
the PGE. We will not actually perform any calculations in this section. Rather 
we will utilize the description of the hxed endpoint solution in the above section 
in conjunction with well known results and certain ideas of Pyne and Birkinshaw 
(1993) to arrive at a plausible understanding. 

Briefly, we set up the mathematics of our lens system. For a Minkowski space 
background with metric = diag(—1,1,1,1), equation (53) takes the form 


(As) = (Ai) + (As - Ai) (Ai) + (As - A) /^^^^(A) dX (56) 

JXi 


An appropriate perturbation for our lens can be determined in the well known 
weak held approximation (see, e.g. Weinberg 1972, Chapter 10) with the result 
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hfjtv = diag(—2(/), —2(j), —2(f), —2(f>) with (f> the Newtonian potential of the lens on 
the spatial hypersnrfaces. It is easiest to snppose that both the emitter and the 
observer are snfficiently far from the lens that its gravitational held npon them may 
be neglected. For simplicity we conld also assnme that both the emitter and the 
observer have fonr-velocities given by = (1, 0, 0, 0). 

Consider the one parameter family of nnll geodesics of onr Minkowski space 
backgronnd given by 




( ^ \ 

2L — 2L cos 9 — \ cos 9 

2Lsin9 + A sin 6* 

0 ; 


(67) 


At this point, there is nothing mathematically to stop ns from constrncting one 
nnll geodesic in the pertnrbed spacetime from each of these backgronnd geodesics, 
snbject to the bonndary conditions x^^^^{—2L) = 0 and In 

addition, we conld hnd by forcing onr constrncted geodesics to be nnll in 

These bonndary conditions wonld ensnre that each of the constrncted 
nnll geodesics intercepts both q and the worldline of the observer. Finally, for each 
pertnrbed geodesic we conld determine the angle it dehned at the observer with the 
image of the lens. 

Of conrse, onr ability to constrnct so many snch geodesics is an enormons warn¬ 
ing sign: we know very well that at most two can be good approximations to the 
exact solntion, one passing to each side of the lens. How then can we single ont the 
two good approximate geodesics from the many bad ones? In Fignre 2, we show 
three of the backgronnd paths we are considering, projected into the xy-plane. One, 
travels diametrically away from the lens, one, travels extremely close to 

the Schwarschild radins of the lens, Rg, and one is that backgronnd path which gen¬ 
erates the best approximate pertnrbed path. The best approximation backgronnd 
path, dehnes an angle with the x-axis at the emitter eqnal to 9e- If we take 

serionsly the cinematic description of onr method offered in section 5, we shonld 
demand that any backgronnd path be snitable for an initial valne calculation with 
initial value specihed at the observer. This would immediately rule out the use 
of the path which passes very close to the lens (Pyne and Birkinshaw 1993). The 
other extreme path, however, is perfectly suitable for such a calculation so that this 
condition is not sufficiently stringent. 

We could proceed by constructing the approximate paths associated to each 
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of the unperturbed paths of Figure 2. Instead, we will use the analysis of section 
5, above, to obtain important qualitative information about the solutions. For a 
given background geodesic, that analysis instructs us, hrst, to consider an initial 
value type solution, corresponding to a photon emitted backwards in time from 
the observer in the direction of the background geodesic. We know from Pyne 
and Birkinshaw (1993) that the initial value calculation returns a deflection at the 
lens equal to La, where a is the usual lens deflection angle appropriate to the 
background path used, a = 4m/Lsm9. The method then would correct for the 
distance between and the origin, essentially by (vectorially) subtracting 

this distance from the deflection computed at the lens plane. The result, considered 
as the net linear deflection undergone by the photon, would then determine the 
position of the emitter on the observer’s sky by applying to it the inverse Jacobi 
operator, basically 1/2L. The paths constructed in this manner are shown in Figure 
3. 

Comparison of Figures 2 and 3, however, reveals an important distinction be¬ 
tween the “correct” solution, and the two extreme solutions, and Xq. The 
two extreme solutions sample a lens potential totally unlike that felt by the back¬ 
ground paths they were constructed from. In contrast, the solution marked x^ feels 
essentially the same lens potential as that felt by the background geodesic it was 
constructed from, x^^\ We could, for instance, conjecture that the two paths we 
seek are those generated by the two background paths which minimize 


I \(j)\dX- I \(j)\dX 

a;(0)M 


(58) 


where x^ is the perturbed path associated to Of course, (58) is mostly heuris¬ 

tic. We are not proposing this expression as a general error functional, assigning 
some signihcant real number measure of the error involved in using any particular 
background geodesic. It does happen to suffice here, however. 

* Nevertheless, such a functional should not be too difficult to hnd. The error in 
the perturbative geodesic expansion is a result of error in the approximate solution 
to Einstein’s equation itself, which we do not consider, error from linearization of 
the Christoffel sysmbols and error from the truncation of the Taylor expansions used 
to express these quantities along the background geodesic (Pyne and Birkinshaw 
1993). Simple matrix methods and the well known remainder term for Taylor’s 
theorem can furnish crude bounds on these last two sources of error. 
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It is important to recognize that the above is exactly the sort of reasoning which 
must be applied in making rough consistency judgements in an initial value calcula¬ 
tion. In fact, there is an analogous background path freedom in those calculations 
as well: by changing the initial data, an inhnite number of background paths can 
be made to generate an inhnite number of distinct perturbed geodesics all of which 
pass though a given point with the same tangent vector. Of course, there is a unique 
actual geodesic which passes through this point with the given tangent. It is nearly 
always assumed in perturbative calculations that the specihc perturbed geodesic 
under investigation is, if not the best, at least an adequate approximation to the 
actual geodesic. In fact, we often have no more reason to expect this in the common 
initial value cases as we do for the hxed endpoint scenarios. It is a common hope, 
however, that forewarned is largely forearmed. 

Suppose, now, that we have applied some error minimization and gained the 
usual lens solution. Would we have gained the usual lens time-delay formula? The 
answer is no, because the usual time-delay formula (see, e.g. Schneider, Ehlers, and 
Falco 1993, Chapter 4) contains a second order term, the geometric delay term, 
which is proportional to the square of the lens angle, itself a hrst order quantity. 
In fact, as we have already noted, the lens perturbation we are considering, when 
used in (10), immediately produces the Shapiro delay, but contributes nothing else. 
From a perturbation-theoretic point of view, the inclusion of the geometric delay 
term in the standard treatments is quite ad hoc. In Seljak (1994), for instance, the 
spatial components of the lensed photon path are solved for to linear order in the 
perturbing potential. This projected path is then “lifted” into the time domain by 
imposing ds‘^ = 0 at second order. We do not intend this remark as a critique on the 
usual treatments: it is possible to perform the entire calculation at second order and 
show that only the usual terms are important. We bring this issue up in order to 
emphasize that there are astrophysical instances when a priori second order terms 
contribute numerically as importantly as hrst order terms. A nice example of this 
is the recent work of Frieman, Harari, and Surpi (1995). 


7. Summary 

We have presented an equation (29) for the time delay measured by two antennae 
observing a distant source through an adiabatic background of gravitational waves 
on Finstein-de Sitter spacetime. We have used the equation to determine the pattern 
of source motions on the sky induced by a background of gravitational radiation. We 
have also shown how these results may be extended to the curved FRW spacetimes. 
Our results are immediately applicable to situations, such as that considered here, 
for which the standard formulae of Kristian and Sachs (1965) do not apply, and 
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so represent an important contribntion to the theory of observations in pertnrbed 
spacetimes. 
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Figure 1. The geometry of Section 3. In the pertnrbed spacetime, light from a 
sonrce at q travels along the nnll geodesic to an antenna at w'^. The antenna is 
separated from the origin of spatial co-ordinates by the spatial vector ta is the 
line of time translates of vj'a- The nnll geodesic of the backgronnd Xa^^ joins the 
sonrce to a point wa G ta- 

Figure 2. The background paths of Section 6, projected into the xy-plane. Rs is 
the Schwarzschild radius of the lensing mass. The arrows point in the direction of 
photon travel. Be is the Einstein angle of the lens. 

Figure 3. The solution paths for the hxed-endpoint problem along the background 
paths of Figure 2. Rs is the Schwarzschild radius of the lensing mass. The arrows 
point in the direction of photon travel. 9e is the Einstein angle of the lens. 
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